A Novel Theory for High Temperature Superconductors considering Inhomogeneous 

Charge Distributions 
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We present a percolation theory for the high- Tc oxides pseudogap and Tc dependence on the hole 
level. The doping dependent inhomogeneous charge structure is modeled by a distribution which 
may represent the stripe morphology and yield a spatial distribution of local Tc{r). The temperature 
onset of spatial dependent superconducting gap is identified with the vanishing of the pseudogap 
temperature T*. The transition to a superconducting state corresponds to the percolation threshold 
among regions of different Tc. As a paradigm we use a Hubbard Hamiltonian with a mean field 
approximation to yield a doping and temperature dependent superconducting d-wave gap. We show 
here that this new approach reproduces the phase diagram, explains and gives new insights on 
several experimental features of high- Tc oxides. 



Pacs Numbers:74.72.-h, 74.20.-z, 74.80.-g, 71.38.-|-i 



I. INTRODUCTION 

High- Tc oxides has been discovered fifteen years agcfl 
but many of their important properties remains not well 
understood. Among these, the pseudogap phenomenon, 
that is, a discrete structure of the energy spectrum above 
Tc, identified by several different experimentstl, has its 
nature not yet been clarified. Such open problem has 
attracted a lot of experimental and theoretical effort be- 
cause it is general belief that its solution is related to 
the understanding of the superconducting fundamental 
interaction. 

The evidence of such energy gap above the supercoiu 
ducting phase is clearly demonstrated by tuaBfilingJ'cl 
and angle-resolved photoemission spectroscopyOfl exper- 
iments. In the resistivity measurements its preaejape is 
seen by a decrease in the linear behavior below T*[]'EI and 
in the specific heat as a suppression in the linear coeffi- 
cient 7 of the temperatureO. 

There is also mounting experimental evidence that the 
hole doped inhomogeneity into the Cu02 planes, com- 
mon to all cuprates, is directly related to the pseudo- 
gap phenomenon. In a given family, the underdoped 
compounds near the doping onset of superconductivity 
have the more inhomogeneous charge distributions and 
the larger T* . As the doping level increases, thej-saiii.- 
ples become more homogeneous while T* decreasegl3il3. 
For overdoped compounds T* disappears or becomes 
equal to Tc- The inhomogeneities were long supposed 
to be importanL. in the studies of high temperature 
superconductorstj pbpt only after the discovered of the 
spin-charge stripealj, they have become matter of sys- 
tematical studies. In the spin-charge stripes scenario, 
regions of the plane are heavily doped (the stripes) and 
other regions are underdoped and fill the space between 
the charge-rich stripes. 

Recently, magnetic excitations as the vortex-like 



Nernst effect have been reported above and a lo- 
cal Meissner state, which usually appears only in the su- 
perconducting phase, has been seen as a precursor to 
superconductivityt3. Such inhomogeneous diamagnetic 
domains develop near T* and grow continuously as T 
decreases towards T^. Near Tc, the domaiiis appear to 
percolate, according Fig. 3 from Iguchi et aO. 

Based on all these experimental findings, about the 
pseudogap phenomenon, local charge inhomogeneities 
and a non percolative local Meissner state between T* 
and Tc, we propose a new scenario to explain the high 
Tc superconductors phenomenology: due to the doping 
dependent charge inhomogeneities in a given compound 
with average charge density pm, there is a distribution 
of local clusters with spatial dependent charge density 
p(r), each with its superconducting transition temper- 
ature Tc(r). T* is the maximum of all Tc(r). As the 
temperature falls below T* , some clusters become super- 
conducting, but they are surrounded by metallic and/or 
antiferromagnetic insulating domains and, consequently, 
the whole system is not a superconductor. The num- 
ber of superconducting clusters increases as the tempera- 
ture decreases, so the superconducting regions grow and, 
eventually, at a temperature Tc, they percolate and be- 
come able to hold a macroscopic dissipationless current. 
Exactly as the Meissner state domains shown in Fig. 3 of 
Iguchi et aO. Similar ideas were discussed by Ovchin- 
nikov at alll3. They were concerned mainly with the 
microscopic mechanism which leads to a distribution of 
Tc(r) and its effect on the density of state. 

In order to show that these ideas are able to make 
quantitative predictions and reproduce the high- Tc ox- 
ides phase diagram, we have performed calculations on 
a Hubbard model and a gap equation is obtained within 
mean field approach. 
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II. THE CHARGE DISTRIBUTION 

The consequence of the microscopic charge inhomo- 
geneities distribution in the Cu02 planes, possibly in a 
striped configuration, is the existence of two phases which 
are spontaneously created in the Cu — O2 planes; regions 
which are heavily doped or hole-rich form the stripes and 
others regions which are hole-poor and are created be- 
tween the charge-rich stripes. The exactly form of these 
charge distributipasjs not known and it is presently mat- 
ter of researchlljO'tS. We have chosen a distribution ca- 
pable to reproduce the experimental observations, and 
for this purpose, we use a combination of a Poisson and 
a Gaussian distribution. For a given compound with an 
average charge density the hole distribution is func- 
tion of the local hole density p, P[p', Pm) divided in two 
branches. The low density branch represents the hole- 
poor or non-conducting regions and the high density one 
represents the hole-rich or metallic regions. As concerns 
the superconductivity, only the properties of the hole rich 
branch are important since the current flows through the 
metallic region. Such distribution may be given by: 

P{p) = ±{p - pc)exp{-{p - p,)V2(a±)2)/ 

(a2)(2-exp(-(0.05)V2(a±)2)) (1) 

The plus sign is for the hole-rich {pc ~ pm) for pm < p, 
the minus to the hole-poor branch {pc — 0.05) for p < 
0.05 and P{p) = for 0.05 < p < Pm- The half-width 
a is related with the degree of inhomogeneities and de- 
creases with . thp ho le density in order to represent current 
observationsHEla. 







Poisson-Gaussian Distribution 




=0.185 Pp=0.542 


P„ 


=0.320 Pp=0.465 


P„ 


=0.520 pp=0.528 (Gaussian) 




/> 

/ \ 

i \ 
i \ 

/'\ ' \ 




1 \ I \ 
\ /' \ 







^ , ,~-r~^ I 

0.2 0.4 '\ if 0.6 0.8 



P 

FIG. 1. Model charge distribution for the inhomogeneities 
or stripe two phase regions. The low density insulating (an- 
tiferromagnetic) branch. The high density hole-rich region 
starts at the compound average density pm and pp, indicat- 
ing by the arrows, is the density where percolation can occur. 

An example of the distribution is shown in fig.l. For 
illustration purpose, we show the results for compounds 



with Pm = 0.185 and (t+ = 0.05 and with = 0.32 
and (7+ = 0.038. Above pm = 0.25 the charge distri- 
bution becomes a simple Gaussian centered at pm with 
a < 0.02, which reflects the non-existence of the stripes 
phases and the observed homogeneous charge distribu- 
tion for the overdoped compounds. 

The values of a for a given sample are chosen in or- 
der that percolation in the hole-rich branch occurs ex- 
actly at a given density pp. Thus Tc{pp) is the maxi- 
mum temperature which the system can percolate and 
which we identify as equal to Tc{pm)- Although we used 
a set of parameters to compare with the experimental 
phase diagram of Bi2212, the main physical aspects can 
be modeled by others distributions with the similar re- 
sults. According to percolation theory, percolation oc- 
curs in a-^quare lattice when 59% of the sites or bonds 
are fiUedta. Thus, we find the density where the hole-rich 
branch percolates integrating J P{p)dp from pm till the 
integral reaches the value of 0.59, where we define pp. Be- 
low Tc{pm) the system percolates and, consequently, it is 
able to hold a dissipationless supercurrent. To estimate 
Tc{pm) we need to calculate T* as function of p. 



III. THE PHASE DIAGRAM 

To develop the dynamics of the hole- type carriers in the 
Cu-0 planes, we adopt a two dimension extended Hub- 
bard Hamiltonian in a square lattice of lattice parameter 
a 

H = - J2 ^ij4aC^J'y + ^ 1] "*T"U 

+ ^ij4a(^]a'(^j'y'(^^'y (2) 

where tij is the hopping integral between sites i and j; 
U is the Coulomb on-site correlated repulsion and Vij 
is the a phenomenological attractive interaction between 
nearest-neighbor sites i and j which will will argue later 
about its possible origin. 

Using a BCS-type mean-field approximation to develop 
Eq.(||) in the momentum space, one obtainSj-the. self- 
consistent gap equation, at finite temperaturesc2rE3 

k' 

with i?k = ^ej^ + Aj^, which contains the attractive 
potential Vkk' in the extended Hubbard Hamiltonian of 
Eq.2. For a d-wave order parameter, U is summed out 
of the gap equation and the amplitude of the attractive 
potential V is the only -unknown variable and become 
an adjustable parametertSo. £k is a dispersion relation. 
The hole density and the gap equation, Mft-solved self- 
consistently for a d-wave order parameterESc^. 
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In Fig. 2 we plot the temperatures of vanishing gap 
from Eq.2 which we take as T*, as function of p. Here we 
use V — — O.lSOeV which reproduces well the experimen- 
tal values of T* X p for the Bi2212 systemlZl. The param- 
eters are -adthin 10% to those taken from ARPES mea- 
surementsEll with a nearest neighbor hopping t = 0.1 2cV 
and further hopping parameters up to fifth neighbor. The 
particular form of the T* curve depends on the values of 
these hoppings. Notice that the density of holes Jaas a 
factor of 2 with respe(jt|to the values given by Ref.Ll, but 
in agreement with rcf.tJ as it is more appropriate for the 
Bi2212 systemN. 
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FIG. 2. The phase diagram taking T*(p) as the onset of 
vanishing gap and Tc{pm) as the percolating threshold. Trtie 
experimental points and the symbols are taking from Ref.Q 

Thus the T* and Tc shown in Fig. 2, for a sample with 
average charge density pm, are obtained in the follow 
way: the hole-rich or metallic branch of the distribution 
describes the regions with hole charge densities p > Pm- 
These charge fluctuations yield clusters with local super- 
conducting temperature Tc{p), as p varies in the sample 
and one may write T'c(r), where r represents any position 
inside the compound. For the metaUic regions, Tc{p) is 
a decreasing function of p and the maximum gap tem- 
perature occurs for T*{pm) = T* . The different metal- 
lic regions in this sample have Tc{p) < T*. For tem- 
peratures below r*, some superconducting clusters are 
formed, like small superconducting islands embedding in 
a metallic and insulating medium. Thus, as the temper- 
ature decreases, more clusters become superconducting, 
and eventually the superconducting regions percolates at 
Tc, that is, a superconducting current can flow for tem- 
peratures T <Tc. 



IV. DISCUSSION 

The fact that T* decreases continuously with p in the 
superconducting region {pm > 0.15), as seen in Fig. 2, 



is very suggestive and in agreement with the early ideas 
regarding a phonon mediated superconducting interac- 
tions: Materials whose vibrating atoms interact strongly 
with the electrons, and are poor metals, should become 
superconductors at higher temperatures than those goiad 
metals, whose atoms interact weakly with electronso. 
For cuprates, as the doping level of the samples increases, 
it is well known that the compounds change from very 
poor metals in the normal phase to very good metals 
with typical Fermi liquid behavior for overdoped sam- 
ples. Since T* x p is a decreasing curve, for any cuprate 
family, and assuming that T* is the onset of supercon- 
ducting gap, such curve may be a strong indication of 
the phononic superconducting interaction. 

There are several observations and measurements that 
can be well explained within the percolating approach, 
we will discuss here only a few examples. 

1-Harris et alii, through ARPES measurements, have 
reported the anomalous behavior of Ao(pm) which de- 
creases steadily with the doping pm although Tc increases 
by a factor of 2 for their underdoped samples. In the over- 
doped region, since Tc also decreases, the behavior is the 
expected conventional proportionality. It is well known 
that superconductors have a constant value for the ratio 
2Ao/A:bTc, being 3.75 for usual isotropic order parameter 
and 4.18 for solutions. 
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FIG. 3. The zero rtamperature gap for 9 samples as mea- 
sured by Harris et alEII and our calculations. 

At low temperature, since the superconducting region 
percolates through different regions, each with a given 
Ao(r), tunneling and ARPES experiments detect the 
largest gap present in the compound. Consequently, 
Ao(pm) must be correlated with the onset of vanishing 
gap T*{pm) which is the largest superconducting temper- 
ature in the sample, and not with Tc{pm)- As we show in 
Fig. 3, correlating the values plotted in Fig. 2 for T*{pm) 
with Ao(pm), we ace, able to give a reasonable fit for the 
data of Harris et alcZl on Dy — BSCCO and explains the 
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different energy scales pointed out by Harris et al and 
several others authors. 

2- The resistivity measurement is also one of the tools 
to detect the pseudogap. The underdoped and optimum 
doped high- Tc oxides have a linear behavior for the re- 
sistivity in the normal phase up to very high tempera- 
ture. However, at T* there is a deviation from the lin- 
ear behavior and.l|he resistivity falls faster with decreas- 
ing temperatureHiH. This behavior can be understood by 
our model, with the increasing of superconducting clus- 
ter numbers and size, as the temperatures decreases be- 
low T*. Each superconducting cluster produces a short 
circuit which decreases the resistivity below the linear 
behavior between T* and Tc- 

3- Recently measurements of magnetic domains above 
Tc has been interpreted as a diamagnetic precursor to 
the Meissner state, produced by peiiormed pairs in un- 
derdoped La2-xSr.j;CuOi thin filmsEj. The existence of 
superconducting cluster between T* and Tc easily ex- 
plains the appearance of local diamagnetic or Meissner 
domains, and, if there is a temperature gradient in the 
sample, the local flux flows and produces the dynamic 
flux flow stateEJ. 

4- Another important consequence is that the pair- 
ing mechanism must be investigated by experiments per- 
formed mainly at T*. A such experiment was accom- 
plished by Rubio Temprano et alcJ, which measured a 
large isotope effect associated with T* and an almost neg- 
ligible isotopic effect associated with Tc in the slightly un- 
derdoped HoBa2Cu^Os compound. The results strongly 
support the fact that electron-phonon induced effects 
are present in the superconducting mechanism associated 
with T* and with the percolation approach to Tc- 



V. CONCLUSIONS 

We have demonstrated that the percolating approach 
for an inhomogencous charge distribution on the Cu02 
planes provides new physical explanations for many ex- 
periments performed on high- Tc oxides and quantitative 
results for their phase diagram. Contrary to some cur- 
rent trends, which Tc is regarded as a phase coherence 
temperature and the existence of a gap phase without 
coherence between Tc and T*, in our approach, Tc is a 
percolating temperature for different regions which, due 
to the inhomogeneities, possess different local supercon- 
ducting transition T*(r). Similarly, instead of having a 
superconducting gap Age and an excitation gap A associ- 
ated with T* , we have a distribution of locally dependent 
Asc{r). 

The method described in this work can be applied 
in any cuprate and yields also several new implications 
which will be discussed elsewhere, but one of the most 
interesting is that one can search for materials with very 
large Tc's if a better control of the local doping level is 
achieved. 
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